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Results for fs, /ss, fo, fog, and their ratios are presented. High statistics quenched runs at /? = 5.7, 5.85, 6.0, 
and 6.3, plus a run still in progress at /? = 6.52 make possible a preliminary extrapolation to the continuum. The 
data allows good control of all systematic errors except for quenching, although not all of the error estimates have 
been finalized. Results from configurations which include effects of dynamical quarks show a significant deviation 
from the quenched results and make possible a crude estimate of the quenching error. 



For the past two years, we have been computing 
heavy-light decay constants with Wilson fermions 
in the quenched approximation, and have more 
recently begun to use Np = 2 dynamical stag- 
gered fermion background configurations to test 
directly the effects of quenching. Computations 
on the largest lattices have been performed on the 
512-node and 1024-node Intel Paragon comput- 
ers at Oak Ridge National Laboratory; Paragons 
at Indiana University and at the San Diego Su- 
percomputer Center have been used for smaller 
lattices. Basic features of the calculation are de- 
scribed in Ref. [1]. Table 1 shows the parameters 
of the lattices used. 

Quenched runs A, E, C, D, and H are on lattices 
with very nearly the same physical volume. Finite 
volume errors in the quenched approximation are 
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determined by comparing the results of runs A 
and B. More work is needed to estimate the finite 
volume errors for the Np = 2 dynamical fermion 
runs. 

In most cases, we compute the pseudoscalar 
decay constant for static-light, as well as heavy- 
light, mesons. However, as explained in [1], the 
current procedure does not produce usable static- 
light results on lattices B, E, I and J. Efforts to 
remedy this situation are in progress. 

For heavy-light mesons we use the Kronfeld- 
Mackenzie norm — 6k) and adjust the mea- 
sured meson pole mass upward by the difference 
between the heavy quark pole mass ("mi") and 
the heavy quark dynamical mass ("012") as cal- 
culated in the tadpole-improved tree approxima- 
tion [2]. Since we do not make the additional 
changes in action and operators called for by the 
heavy Wilson quark program [3], some 0(ma) er- 
rors will still be present. Further 0(a) errors are 



2 



Table 1 

Lattice parameters. Runs F, G, I, J, and K use 
variable-mass Wilson valence quarks and two fla- 
vors of flxed-mass staggered dynamical fermions; 
all other runs are quenched with Wilson valence 
quarks. The numbers in parenthesis are the 
planned flnal number of conflgurations for runs 
still in progress. F lattices were provided by the 
Columbia group; G, by HEMCGC 



name 




size 


# conflgs. 


A 


5.7 


8^ >< 


48 


200 


B 


5.7 


16^ 


>< 48 


100 


E 


5.85 


12^ 


>< 48 


100 


C 


6.0 


16^ 


>< 48 


100 


D 


6.3 


243 


>< 80 


100 


H 


6.52 


32^ >< 


100 


18 (50) 


F 


5.7 


16^ 


>< 32 


49 




m = 0.01 








G 


5.6 


16^ 


>< 32 


95 (200) 




m = 0.01 








I 


5.445 


16^ 


>< 24 


118 




m = 0.025 








J 


5.47 


16^ 


>< 24 


128 




m = 0.05 








K 


5.415 


8^ >< 


24 


78 




m = 0.0125 









introduced by the Wilson light quark. An extrap- 
olation to a = is thus crucial. 

Since we only have results for degenerate light 
quarks, we determine the strange quark hopping 
parameter Kg by adjusting the pseudoscalar mass 
to ^2m|f — , the lowest order chiral pertur- 
bation theory value. 

A plot of fp\/Mp vs. 1/Mp is shown in flg. 1 
for lattice D. The flts are to the form co + ci/Mp + 
C'j/Mp. Because the systematic errors are ex- 
pected to be larger for the heavier propagating 
quarks, we use flts like that shown to the "lighter 
heavies" (with meson masses lighter than the D) 
on each lattice to give our central values for de- 
cay constants. Fits to the "heavier heavies" (with 
meson masses from the D region up to 4 GeV) 
give one estimate of the systematic error caused 
by large lattice masses. 

Figure 2 shows fp plotted versus lattice spac- 
ing, as determined from Extrapolation to 
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Figure 1. fp{Mp)i vs. 1/Mp for lattice D. The 
lower line is a covariant flt (x^/d.o.f. = 1.1) to the 
diamonds ("lighter heavies" -|- static); the slightly 
higher line is a covariant flt (x^/d.o.f. = 2.7) to 
the squares ("heavier heavies" -|- static). The 
light quark is extrapolated to the physical mass 
(niu + md)/2. The scale is set by 

a = of the quenched data by a linear flt seems 
justifled; the flt has x^/d.o.f. = 1.2 and gives 
fp = 151(5)MeV. We are unwilling at this point 
to attempt an extrapolation of the dynamical 
fermion data. The present data from the three 
such runs at the largest values of a (runs I, J and 
K) is suspect since the lattices are short in the 
time direction (Ni = 24) and the plateaus we do 
have are rather poor. 

Work is in progress on longer versions of the 
lattices at large a and on lattices at intermediate 
values of a. We hope that an extrapolation of the 
dynamical fermion data to the continuum will be 
possible when that work is completed. 

The systematic errors are estimated by varying 
the parameters and flts used in the computation. 
Table 2 shows central values and error estimates 
for several decay constants and ratios from vari- 
ous versions of the analysis. Some errors are es- 
timated in more than one way, and not all the 
different errors are in principle independent. The 
versions are as follows: 

0) Central values: Fits to "lighter heavies" with 
scale and linear extrapolation in a. 
Boosted couplings approximately given by 
glfP = (jf^/(plaquette) but with no tadpole 
improvement of the zero-mass axial renor- 
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Figure 2. fs vs. a. Octagons are quenched lat- 
tices (runs A, B, E, C, D, H); crosses are Np = 2 
dynamical fermion lattices (runs F, G, I, J, K). 
A linear fit to the quenched points only is shown; 
the extrapolated value at a = is indicated by 
the square. The slightly lower quenched result at 
a K. 0.9GeV~^ comes from run B; the higher one, 
from run A. The scale is set by = 132 MeV 
throughout. 

malization constant Za- Fields normalized 
by \/l — 6k; meson mass shifted upward by 
1112 — rrii ; Kg fixed from the pseudoscalar 
mesons as described above. 

1) Statistical errors in version 0. Experience 

shows that systematic errors due to choice 
of fitting ranges in t of the correlators are 
comparable to the statistical errors. 

2) Errors due to interpolation in l/M at fixed a; 

estimated by replacing correlated l/M fits 
in version by uncorrelated fits. 

3) (^((ma)^) errors: Replace "lighter heavies" in 

version by "heavier heavies." (0(am) er- 
rors are presumably corrected for by the lin- 
ear extrapolation in a.) 

4) 0((ma)^) errors: Instead of ru'j — mi, shift 

the meson mass by nis — nii where nis is 
the "hyperfine splitting" mass (l/nis is the 
coefficient of cr • 5) [2]. 

5) Finite volume errors: Take percent difference 



Table 2 

Central values and errors from various versions 
of the analysis for fs, fsj/s, fo,, and foj Id- 
See text for description of the versions. Decay 
constants and errors are in MeV. 



version 


fB 






Id J Id 


0) 


151 


1.11 


198 


1.09 


1) 


5 


0.02 


5 


0.01 


2) 


10 


0.00 


3 


0.01 


3) 


1 


0.02 


2 


0.00 


4) 


5 


0.00 


3 


0.01 


5) 


8 


0.01 


4 


0.01 


6) 


4 


0.01 


6 


0.03 


7) 


5 


0.01 


3 


0.00 


8) 




0.07 


7 


0.05 


9) 


11 


0.08 


1 


0.04 


10) 


26 


0.00 


14 


0.01 


11) 




0.03 


19 


0.02 



between runs A and B and apply that per- 
centage to final, extrapolated values. 

6) 0{a^) errors: Replace linear extrapolation in 

a by linear plus quadratic extrapolation. 

7) Weak coupling perturbation theory error: If 

boosted couplings in version are expressed 
as gviq*) [4] then q* ^ 4.3/a. Change q* 
to 2.2/a (still with no tadpole improvement 
of Za) and compare. 

8) Error from fixing Kg'. Determine Kg by adjust- 

ing mass of vector ss state ((f)) to physical 
value and compare with version 0. 

9) Error from fixing the scale: Use nip instead of 

to set the scale. 

10) Quenching: Compare result from run G to 

the quenched results (a la version 0) in- 
terpolated to the same lattice spacing, and 
take same percentage error on final extrap- 
olated result. 

11) Quenching: Same as version 10, but also fix 

Kg from the (f) everywhere, as in version 8. 
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The systematic error within the quenched ap- 
proximation is now determined by adding in 
quadrature the errors determined by versions 1 
(used as estimate of errors in t fits) + 2-1-5-1-6 
+ 7 + the larger of 3 or 4. Note that if version 
6 were guaranteed to be a good estimator of the 
error in the extrapolation to the continuum, ver- 
sions 3 or 4 would be superfiuous. However, since 
the error in these error estimates is probably com- 
parable to the errors themselves, we take a more 
conservative approach and include the larger of 
versions 3 or 4 in addition to version 6. 

Once the quenched approximation computa- 
tion has been reliably extrapolated to the con- 
tinuum, any dependence of the results on how 
the scale is fixed (version 9) or how Kg is fixed 
(version 8) must be considered as an error due to 
quenching itself. Assuming our extrapolation is 
reliable, we therefore estimate the quenching er- 
ror by taking the largest of the errors determined 
by versions 8, 9, 10, or 11. 

It is clear that some of the estimates given in 
Table 2 are rather crude; all should be taken as 
preliminary. First of all, the use of perturba- 
tion theory should be rationalized by using tad- 
pole improvement throughout and then compar- 
ing several reasonable choices for the scale of the 
coupling constant. Further, the effect of chang- 
ing the fitting ranges in t should be determined 
directly. Errors in the interpolation in l/M are 
more naturally determined by changing the fit- 
ting function rather than the type of fit. C. All- 
ton has suggested [7] estimating the large am er- 
rors by comparing to an "old-style" analysis with 
\/2k normalization and no mass shifts. He has 
further advocated estimating the error in the con- 
tinuum extrapolation by fitting the three points 
with smallest a to a constant. These approaches 
will be tried. 

Finally, the quenching error can be more re- 
liably estimated when the Np = 2 dynamical 
fermion results can be extrapolated to the contin- 
uum. Note from fig. 2 that the error determined 
thereby may be considerably larger than the cur- 
rent estimate at fixed lattice spacing. Of course, 
at that point one may wish to take as central val- 
ues the Np = 2 results. However, there would 
still be a quenching error since one is (at least at 



present) only extrapolating the valence and not 
the dynamical quarks to the chiral limit. (This 
is a "partially quenched theory" [5].) The real 
world, moreover, has Np = 3. In fact the Np = 2 
simulations, even after extrapolation to the con- 
tinuum, may be just as far from the real world as 
those of the quenched approximation [6]. 
The (still preliminary) results are: 

fB = 151(5)(16)(26) fo = 182(3)(9)(22) 
/s, = 169(7)(14)(29) fo, = 198(5)(10)(19) 

^ = 1.11(2)(4)(8) ^ = 1.09(1)(4)(5) 

JD 

where the first error is statistical; the second, 
the systematic error within the quenched approx- 
imation; the third, the estimate of the quench- 
ing error. Decay constants are in MeV. While 
the estimate of the quenching error has large 
uncertainties at this stage, the error within the 
quenched approximation will, we think, be deter- 
mined rather well once the afore-mentioned re- 
finements in the analysis are completed. 
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